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Abstract 

We calculate the near-field optical spectra of excitons and biexcitons in semiconductor quantum 
dots naturally occurring at interface fluctuations in GaAs-based quantum wells, using a non-local 
description of the response function to a spatially modulated electro-magnetic field. The relative 
intensity of the lowest, far-field forbidden excitonic states is predicted; the spatial extension of the 
ground biexciton state is found in agreement with recently published experiments. 

PACS numbers: 73.21.La,78.67.-n,71.35.-y 
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Optical manipulation of states in quantum dots (QDs) is a key towards the implementa- 
tion of quantum information processing in semiconductors. ^'^'^ Important features of QDs 
are not only that excitations are fully discrete in energy resulting in long coherence times, ^'^ 
but also that they are governed by the Coulomb interaction between electrons and holes 
which gives rise to specific few-particle aggregates as a result of photoexcitation. Based on 
this unique mechanism, all-optical conditional control of exciton and biexciton states were 
proposed and demonstrated for quantum computing schemes in QD structures. ^'^ 

The use of local probes for a detailed understanding of these interacting states and their 
real-space distribution is thus of great interest for their engineering, and might open the way 
to the important goal of direct optical manipulation of individual states in coupled struc- 
tures. The so-called natural QDs — or terraces — , where quantum confinement is induced by 
local monolayer fluctuations in the thickness of a semiconductor quantum well, have sofar 
been excellent laboratories for the study of these novel phenomena because quantization is 
combined with a large oscillator strength associated to the large active volume of the QD.''' 
The huge progress in probe preparation has lead to resolutions which in this class of systems 
have reached the scale of the constituent quantum states. 

In this letter, we investigate theoretically the optical response of excitons and biexcitons 
in natural QDs in near-field experiments. We calculate the spatial maps of the exciton and 
biexciton low-energy states, taking into account the non-local susceptibility of the electro- 
magnetic (EM) field, and show that the spatial extension depends on the correlated state, 
being narrower for biexcitons, in agreement with experimental findings. We also calculate 
the oscillator strengths of the lowest, far-field forbidden excitonic states, and show that they 
can be comparable to the ground state, thus allowing for their observation at the presently 
achievable resolution scale. 

Exciton and biexciton states. For excitonic complexes bound to interface fiuctuations, 
where the confinement considerably exceeds the excitonic Bohr radius, center-of-mass and 
relative motion of the electron-hole pairs can be decoupled. Then, the in-plane part of the 
exciton and biexciton wavef unctions, and respectively, can be approximately written 
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Here, 0^ and 0^ are the exciton and biexciton wavefunctions for an ideal quantum well, 
whereas and <i>^ are the envelope functions accounting for interface fluctuations. We 
denote electrons with Vi and r2, holes with Va and r^, and finally the exciton and biex- 



taken from Kleinman^^ for a 5 nm thick quantum well, i.e., the usual two-dimensional ex- 
citon state oc exp— A;ria, k = 1.94 in units of the Bohr radius Go, and the biexciton state 
'^{ria,f^ib,'>^2a,i^2b)x{f^ab), which cousists of trial wavefunctions if) and x associated to the at- 
tractive electron-hole^^ and repulsive hole-hole^^ interactions, respectively; r^j is the distance 
between particles i and j. 

We consider a prototypical interface-fluctuation confinement of rectangular shape with 
dimensions 100 x 70 nm^, and monolayer fluctuations of a 5 nm thick well.^°'^^ The exciton 
confinement is then given by the well-width dependent electron and hole single-particle 
energies along z convoluted with the probabilities Pi{r — R) of finding an electron or hole 
at distance r from the center-of-mass coordinate J2,^^ with a straightforward generalization 
to biexcitons (r and r denote the phase space for electrons and biexcitons, respectively) 



Figure 1 shows the effective confinement potential for excitons (solid line) and biexcitons 
(dashed line). The insets show that for excitons the electron probability distribution Pe{r) 
extends over the effective Bohr radius k~^ao whereas Phif ) is strongly peaked around 0; in 
contrast, for the biexcitons the weaker Coulomb binding results in a strong delocalization of 
the two-exciton complex. As consequence, the biexciton center-of-mass motion is confined 
within a significantly smaller region. 

The resulting two-dimensional Schrodinger equation for the exciton and biexciton center- 
of-mass wavefunctions ^^{R^) and ^^{Rf,), respectively, are solved through a real-space 



citon center-of-mass coordinates with R^ and Ri,. The specific forms of (^o and (^o are 
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discretization and numerical diagonalization on a grid of typical dimensions 64 x 64, similarly 
to the procedure described in Refs. 16. The computed exciton and biexciton states exhibit 
symmetries reminiscent of the two-dimensional box-like confinement,^^ i.e., an s-like exciton 
groundstatc (fig. 2a), two jo-like excited states of lowest energy with nodes along x (fig. 2b) 
and y (not shown), and two nodes along x for the third excited exciton state (fig. 2c); 
finally, the biexciton groundstate (fig. 2d) indeed shows a much stronger localization than 
the exciton groundstate (fig. 2a). 

Near-field spectra. When the EM-field is modulated on the scale of the relevant quantum 
states, the non-local response of the system must be taken into account. In addition to 
the relaxation of the far-field selection rules due to the different symmetry of the field with 
respect to the quantum states, spatial coherence may give rise to interference effects, so 
that collected spatial maps may be non-trivially related to the localization of the excitonic 
wavef unctions. We compute the near-field spectra analogously to the procedure described in 
Refs. 17,18 for an EM-field distribution ^(i^ip — r) of Gaussian shape centered around the 
tip position i^ip. For the excitons the local absorption spectra at a given exciton energy 
are then given by the square modulus of the convolution of ^ with the exciton wavefunction 
^a;_ 16,18 p^j, ^^iq bicxciton we have to be more specific of how the system is excited. We shall 
assume that the QD is initially populated by the ground state exciton and that the near-field 
tip probes the transition to the biexciton states. This situation approximately corresponds 
to that of Ref. 10,11 with non-resonant excitation in the non-linear power regime. The local 
spectra are then proportional to^^ / dr ^(i?tip — f){A'^h{f)i^e{i')\b) , with x and h denoting 
the exciton and biexciton states of cq. (1), respectively, and ■ip^(^h)[r) is the usual fermionic 
field operator for electrons (holes). Then, the square modulus of 



convoluted with ^ gives the optical near-field spectra corresponding to the transition from 
X to b, with ii{r — R) = f dT 6{R — Rx) (poifi, ra)(t>o{r, ri, r, r^) giving the probability of 
exciting an electron-hole pair at r when an exciton is located at R (see fig. le)^^. 

In the second and third rows of fig. 2 we report our calculated optical near-field spectra 
for spatial resolutions of 25 and 50 nm. It should be noted that the first (fig. 2b) and second 
excited state (not shown) are dipole forbidden, but have large oscillator strengths for both 
resolutions. Note also that, as a result of interference effects, the spatial maps at finite 
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spatial resolutions differ somewhat from the wavefunction maps, particularly for the excited 
states: the apparent localization is weaker and, in (c), the central lobe is very weak for both 
resolutions. Finally, we observe that for the smaller spatial resolution the biexciton ground 
state depicts a stronger degree of localization than the exciton one, in nice agreement with 
the recent experiment of Matsuda et al.^°'^^. 

In order to be more quantitative and show if near-field experiments may distinguish the 
dipole forbidden transitions, we have calculated the total absorbed power (i.e., the incoher- 
ently summed intensity of the maps exemplified in fig. 2) as a function of the spatial resolu- 
tion. Figure 3 shows that the intensity of the lowest dipole-forbidden state is a substantial 
fraction of the ground state intensity, and larger than the next dipole allowed excited states 
up to resolutions comparable with the QD linear dimensions. Indeed, preliminary results 
indicate that p-like structures have been seen in the sample of Ref. 9.^° 

We thank K. Matsuda for helpful discussions. This work has been supported in part 
by the EU under the TMR network "Exciting", the Austrian science fund FWF under 
project P15752-N08, the Itahan Minister for University and Research under project FIRB 
Quantum phases of ultra-low electron density semiconductor hetero structures, and by INFM 
I.T. Calcolo Parallelo (2003). 
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With the material parameters of Ref. 18 we find energies of —8.686, —8.442, —8.128, and 
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—8.043 meV for the exciton groundstate and the first three excited exciton states, respectively, 
and —17.264 meV for the biexciton groundstate; exciton (biexciton) energy zero is given by 
the energy of the two-dimensional exciton (biexciton) of the narrow quantum well.-^^ Note that 
we are not attempting a quantitative comparison of the total exciton and biexciton energies 
because of the well-known problems regarding the underestimation of the biexciton binding for 
the trial wavefunction under consideration; for a discussion see, e.g., O. Heller, Ph. Lelong, and 
G. Bastard, Phys. Rev. B 56, 4702 (1997). 
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FIG. 1: Confinement potential along x for the center-of-mass motion of excitons (solid line) and 
biexcitons (dashed line).^^ The insets report the exciton probability distributions Pi{r) for (a) 
electrons and (b) holes, the biexciton probability distributions pi(r) for (c) electrons and (d) holes, 
and (e) /x(r) the probability of creating a second electron-hole pair at distance r from the exciton 
center of mass. The reduced probability at the center of (d) is attributed to the repulsive part x 
of the trial wavefunction.^^ 

FIG. 2: (a-d) Real-space map of the square modulus of the wavefunctions for the exciton (a) 
ground state, (b) first and (c) third excited state, and (d) the biexciton ground state;^^ the dashed 
lines indicate the boundaries of the assumed interface fluctuation, (a'-d') Near-field spectra for a 
spatial resolution of 25 nm and (a"-d") 50 nm, as computed according to Refs. 17,18 and eq. (3). 
The FWHM of the assumed EM-field distribution is indicated in the 2°^^ and 2,^^ row 

FIG. 3: Total absorbed power of the first (black, dashed line), second (gray, dashed line), and third 
(black, solid line) excited exciton state as a function of the spatial resolution of ^. For all resolutions 

the spectra are normalized to the ground state absorption. The first and second excited excitons 
are far-field forbidden, whereas the oscillator strengths of the third excited state is approximately 
one ninth of that of the ground state. 
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